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ABSTRACT 


This  implementation  of  the  physical  optics  approximation  method  for  predicting 
the  radar  cross  section  (RCS)  of  complex  objects  utilizes  the  scientific  computational 
features  of  MATLAB  and  its  Graphical  User  Interface  (GUI)  functions  to  provide  an 
error-free  encoding  of  input  parameters  and  efficient  calculation.  The  POFACETS  GUI 
provides  a  convenient  tool  for  both  students  and  electromagnetic  professionals  to  design 
complex  models  by  representing  its  component  parts  as  triangular  facets  and  offers 
options  for  defining  the  surface  roughness  of  the  scattering  object.  The  GUI  calculates 
the  monostatic  or  bistatic  radar  cross  section  of  the  modeled  object  based  on  certain  radar 
frequencies  and  parameters  of  interest,  creates  a  variety  of  plots  for  the  model  geometry 
and  its  RCS  measurement,  and  provides  additional  functionalities  to  enhance 
visualization  of  the  object  and  verify  correctness  of  data.  Various  help  windows  and 
messages  are  available  to  the  user  to  serve  as  guide  in  the  efficient  use  of  the  GUI.  The 
simplicity  of  the  GUI  design  does  not  require  the  user  to  have  a  detailed  understanding  of 
the  various  steps  involved  in  the  RCS  calculation. 
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I.  INTRODUCTION 


Designing  operational  platforms  and  predicting  their  radar  cross  section  (RCS)  for 
a  range  of  frequencies  provides  a  challenging  problem  in  modem  electronic  warfare.  One 
approach  to  high  frequency  prediction  calculations  is  to  estimate  a  complex  model  with 
an  array  of  simple  shapes,  such  as  triangles  or  flat  plates.  The  RCS  is  obtained  by 
computing  the  scattered  field  of  the  collection  of  these  simple  shapes  to  obtain  the  total 
RCS  of  the  target. 

Among  the  most  common  RCS  prediction  methods  for  any  arbitrary  three- 
dimensional  target  are  the  Method  of  Moments  (MOM),  the  Finite  Difference  Method 
(FDM),  Geometrical  Optics  (GO)  and  Physical  Optics  (PO).  The  MOM  is  a  popular 
technique  that  solves  the  electric  or  magnetic  field  integral  equation  by  reducing  it  to  a 
matrix  problem  with  the  matrix  size  directly  related  to  the  electrical  size  of  the  scattering 
body  in  terms  of  the  incident  wave.  This  method  gives  very  accurate  results  but  is 
computationally  intense.  In  the  FDM,  finite  differences  are  used  to  approximate  the 
differential  operators  in  Maxwell’s  equations  in  either  the  time  or  frequency  domain.  The 
target  must  be  discretized  as  in  the  MOM.  Since  this  method  necessarily  calculates  the 
fields  in  the  region  around  the  target,  computing  the  RCS  of  a  target  with  a  characteristic 
dimension  of  several  orders  of  magnitude  of  the  wavelength  would  entail  considerable 
time.  GO  and  the  Geometric  Theory  of  Diffraction  (GTD)  involve  ray  tracing  methods  to 
analyze  electrically  large  targets  of  arbitrary  shape.  Diffracted  rays  are  allowed  to 
originate  from  the  scattering  of  the  incident  wave  at  edges,  comers,  and  vertices.  Among 
the  problems  associated  with  GO  is  that  the  source  and  the  observation  point  must  satisfy 
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Snell’s  Law  or  the  reflected  field  is  zero.  It  yields  discontinuous  fields  at  shadow 
boundaries  and  zero  fields  in  the  shadow  regions.  GTD  solves  this  problem  by 
accounting  for  the  edge-related  effects.  PO  overcomes  the  problems  associated  with  the 
infinities  of  ray  tracing  for  flat  and  singly  curved  surfaces  by  approximating  the  induced 
surface  fields.  The  induced  currents  are  integrated  over  the  illuminated  portions  of  the 
target  to  obtain  the  scattered  far  field  while  setting  the  current  to  zero  over  the  shadowed 
portions.  It  is  most  accurate  in  the  specular  direction  and  gives  good  estimates  at  points 
far  from  discontinuities  and  shadow  boundaries.  However,  because  PO  abruptly  sets  the 
current  to  zero  at  a  shadow  boundary,  the  computed  field  values  at  wide  angles  and  in  the 
shadow  regions  are  inaccurate.  Also,  surface  waves,  multiple  reflections  and  edge 
diffractions  are  not  included  [Refs.  1, 2  and  3]. 

There  is  a  need  for  an  easy-to-use  tool  that  conveniently  allows  the  modeling  of 
an  arbitrary  three-dimensional  object  with  component  triangular  facets,  visualization  of 
the  geometry,  and  calculation  of  the  radar  cross  section.  The  physical  optics 
approximation  is  an  efficient  method  that  gives  accurate  results  for  most  large  targets,  yet 
is  not  too  computationally  demanding.  To  make  this  tool  appealing  to  the  user  and 
prevent,  if  not  eliminate  unnecessary  errors  in  the  setting  up  of  various  parameters,  the 
design  of  a  graphical  user  interface  is  desirable. 

This  thesis  focuses  on  the  design,  development,  and  implementation  of  a 
Graphical  User  Interface  (GUI)  in  MATLAB  for  a  PO-based  RCS  prediction  code.  The 
program  allows  creation  of  a  model  comprised  of  triangular  facets  with  options  for 
specifying  the  surface  characteristics.  It  then  calculates  the  monostatic  or  bistatic  radar 
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cross  sections  given  certain  computational  parameters  without  requiring  the  user  to  have 
a  detailed  understanding  of  the  various  steps  involved  in  the  calculation. 

The  end  product  of  this  thesis  is  the  POFACETS  GUI  that  implements  a  physical 
optics  approximation  method  that  may  be  used  by  both  students  and  electromagnetic 
professionals  as  a  convenient  tool  to  design  complex  models  and  predict  their  RCS 
performance  based  on  the  radar  frequency  and  other  parameters  of  interest.  The  GUI 
provides  additional  functionalities  to  enhance  visualization  of  the  model  geometry  and 
generates  a  variety  of  plots  to  present  the  RCS  data.  The  GUI  is  based  on  codes 
developed  by  Professor  David  Jenn  and  utilizes  the  scientific  computational  features  of 
MATLAB. 

Chapter  II  gives  a  brief  description  of  the  theory  behind  the  formulation  of  the 
code  by  providing  an  overview  of  radar  fundamentals,  electromagnetic  scattering  and  a 
comparative  discussion  of  the  common  methods  of  high  frequency  RCS  prediction. 

In  Chapter  in  the  flow  or  framework  of  the  code  is  presented  with  the  supporting 
mathematics  involved  in  using  the  physical  optics  method  to  calculate  the  RCS. 

Chapter  IV  describes  the  POFACETS  Graphical  User  Interface.  It  is  intended  to 
provide  a  descriptive  guide  in  making  the  user  acquainted  with  the  operation  of  the  GUI. 
A  step-by-step  illustration  is  provided  for  working  on  the  examples  that  are  discussed. 

Finally,  Chapter  V  summarizes  the  research  and  suggests  further  work  that  could 
be  performed  to  enhance  the  code. 
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II.  THEORETICAL  FORMULATION 


A.  RADAR  BACKGROUND 

The  acronym  radar  was  coined  during  World  War  II  from  the  phrase  radio 
detection  and  ranging.  A  radar  is  a  device  that  transmits  an  electromagnetic  wave  and 
detects  objects  by  virtue  of  the  energy  scattered  from  them  in  the  direction  of  the 
receiver.  Generally,  the  radar  transmitter  and  receiver  can  be  at  different  locations  when 
viewed  from  the  target  and  this  is  called  a  bistatic  radar.  In  most  applications  the 
transmitter  and  receiver  are  co-located  and  often  share  the  same  antenna,  a  configuration 
referred  to  as  a  monostatic  radar.  The  classical  form  of  the  radar  equation  is 


Pr  =  (•£%)(— 
4nR  4 7tR 


^  A  rrr 


where  Pr  is  the  power  received  (in  watts),  Pt  the  output  power  of  the  transmitter,  Gt  the 
gain  of  the  transmitter,  Gr  the  gain  of  the  receiver,  crthe  radar  cross  section  (in  square 


meters),  A  the  wavelength  (in  meters),  L  the  reduction  in  the  received  signal  due  to  losses 
(0  <  L  <  1) ,  and  R  the  range  of  the  target  from  the  radar.  The  first  term  in  parenthesis  is 
the  power  density  at  the  target  (in  watts/m2)  due  to  radiation  from  the  transmitter.  The 
second  term  is  the  incident  power  reflected  or  scattered  back  to  the  receiver.  The  third 
term  is  the  amount  of  the  reflected  power  which  is  captured  by  the  receiving  antenna  due 
to  its  effective  area.  For  the  case  of  a  monostatic  radar  with  the  same  antenna  used  for 


transmitting  and  receiving  (G(is  the  same  as  Gr ),  equation  2.1  simplifies  to 

P&oXh 
r  ~  (4 7lfR* 


(2.2) 
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B.  ELECTROMAGNETIC  SCATTERING 


1.  Radar  Cross  Section 

Radar  cross  section  (RCS)  is  a  measure  of  the  power  scattered  in  a  given  direction 
when  a  target  is  illuminated  by  an  incident  wave,  normalized  to  the  power  density  of  the 
incident  field.  The  purpose  of  the  normalization  is  to  remove  the  effect  of  the  range  and 
arrive  at  an  RCS  computation  that  is  independent  of  the  distance  between  the  target  and 

illumination  source.  In  terms  of  the  incident  and  scattered  electric  field  intensities,  Et 

and  Es ,  the  RCS  is  defined  as 


cr=  iim  4 nR2 

R  — )■ 


(2.3) 


I  *1 

The  unit  of  RCS  most  commonly  used  is  decibels  relative  to  a  square  meter  (dBsm): 

cr,  dBsm  =  101og(<j,m2)  (2.4) 

Radar  cross  section  is  a  function  of  many  factors  which  include  the  target  configuration 
and  its  material  composition,  frequency  or  wavelength,  transmitter  and  receiver 
polarizations,  and  the  target  aspect  (angular  orientation  of  the  target)  relative  to  the  radar. 
Hence,  in  general,  a  could  be  specified  as  apq  ( 6,0)  where  q  and  p  pertain  to  the  incident 
and  received  polarizations,  respectively,  and  where  ( 6, 0)  are  the  polar  spherical  angles  of 
view  as  shown  in  Figure  1. 
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Figure  1.  Coordinate  System  for  Incident  Polarization 


2.  Scattering  Regions 

The  scattering  characteristics  of  a  target  are  dependent  on  the  frequency  of  the 
incident  wave.  Three  regimes  distinctly  differentiate  the  RCS  scattering  of  a  target 
depending  on  the  ratio  of  the  wavelength  X  to  the  body  size  L.  The  first  is  the  low 
frequency  region  or  Rayleigh  region  (A  »  L)  where  the  phase  variation  of  the  incident 
plane  wave  across  the  target  is  small  and  the  shape  of  the  body  is  not  of  importance.  The 
second  is  the  resonance  region  or  Mie  region  (A  —  L)  where  the  phase  variation  of  the 
current  across  the  body  is  significant  and  all  parts  contribute  to  the  scattering  pattern.  The 
third  is  the  high  frequency  region  or  Optical  region  (A  «  L)  where  the  phase  of  the 
current  across  the  body  varies  with  many  cycles  and  the  scattering  process  can  often  be 
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considered  as  the  summation  of  the  returns  from  non-interacting  isolated  scattering 
centers  [Refs.  1, 2  and  4], 

3.  Target  Scattering  Matrix 

The  radar  cross  section  is  a  scalar  number  that  is  a  function  of  the  polarization  of 
the  incident  and  received  waves.  The  polarization  of  the  incident  wave  from  the  radar 
will  not  always  be  in  the  same  orientation  as  the  polarization  of  the  target  coordinate 
system  and  thus,  it  is  necessary  to  decompose  the  polarization  of  the  incident  wave  into 
components  in  the  target  coordinate  system.  Since  the  wave  vector  components  at  large 
distances  are  tangent  to  a  sphere,  two  orthogonal  components  in  terms  of  the  variables  6 
and  <f>  are  sufficient  to  represent  the  incident  field  in  a  spherical  system  centered  at  the 
target.  Thus,  in  general,  the  incident  field  can  be  written  as 

E,  =  Ewd  +  (2.5) 

A  A 

where  6  and  0  are  unit  vectors  in  the  target  coordinate  system.  The  polarization  of  the 
scattered  field  is  not  necessarily  the  same  as  the  incident  field  since  most  complex  targets 
generate  a  cross-polarized  scattering  component  due  to  multiple  diffractions  and 
reflections.  The  following  scattering  matrix  can  be  used  to  specify  the  relationship 
between  the  incident  and  scattered  fields: 

Ese  See^etp  Eig 

E  S  S  E  (2*6) 

s<j>  J  _  00  ^  00  ^i0 

The  Spq  are  the  scattering  parameters,  where  the  first  index  specifies  the 

polarization  of  the  receive  antenna  and  the  second,  the  polarization  of  the  incident  wave. 
The  elements  of  the  scattering  matrix  are  complex  quantities  and  in  terms  of  the  RCS 
[Ref.  1], 
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(2.7) 


S 


pq 


JEH  j<r* 

y[4ftR 2 


=  tan 


'¥W' 

^Re(£jp /£,,)“ 


B.  COMMON  RCS  PREDICTION  METHODS 

Among  the  most  common  numerical  RCS  prediction  methods  for  any  arbitrary 
three-dimensional  target  are  the  Method  of  Moments  (MOM),  the  Finite  Difference 
Method  (FDM),  Geometrical  Optics  (GO)  and  Physical  Optics  (PO).  The  first  two 
provide  rigorous  solutions  but  are  computationally  intense  requiring  considerable 
computer  time  to  calculate  results,  while  the  latter  two  are  approximate  methods  which 
give  accurate  and  fast  results,  but  have  limitations. 

1.  Method  of  Moments 

The  MOM  involves  breaking  the  target  into  a  number  of  surface  or  line 
subdomains.  This  popular  solution  technique  solves  the  electric  or  magnetic  field 
integral  equation  by  reducing  it  to  a  matrix  problem  with  the  matrix  size  directly  related 
to  the  electrical  size  of  the  scattering  body  in  terms  of  the  incident  wavelength.  The 
coefficients  of  a  series  approximation  to  the  current  can  then  be  determined  using  matrix 
algebra.  In  general,  the  subdomains  must  be  small  compared  to  the  wavelength.  This 
method  gives  very  accurate  solutions  and  is  used  in  programs  such  as  the  Numerical 
Electromagnetics  Code  (NEC)  [Ref.  3].  However,  the  calculation  time  is  dictated  by  that 
required  to  solve  the  simultaneous  equations  and  the  method  may  not  be  practical  for 
large  targets  at  high  frequencies  [Refs.  1  and  3]. 
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2. 


Finite  Difference  Methods 


In  this  method,  finite  differences  are  used  to  approximate  the  differential 
operators  in  Maxwell’s  equations,  usually  in  the  time  domain.  The  target  must  be 
discretized  as  in  the  MOM  and  the  elements  have  to  be  small  compared  to  the 
wavelength.  Finite  difference  techniques  do  not  require  the  solution  of  large  dense 
matrices  that  the  MOM  does  because  the  solution  is  stepped  in  time  throughout  the 
scattering  body.  Since  this  method  calculates  the  fields  in  a  computational  grid  around 
the  target,  calculating  the  RCS  of  a  target  with  a  characteristic  dimension  of  several 
orders  of  magnitude  of  the  wavelength  would  entail  considerable  time  to  compute  [Refs. 
1  and  3]. 

3.  Geometrical  Optics 

GO  involves  ray  tracing  methods  to  analyze  electrically  large  targets  of  arbitrary 
shape.  The  rules  for  ray  tracing  in  a  simple  medium  (linear,  homogeneous,  and  isotropic) 
are  similar  to  reflection  and  refraction  in  optics.  In  addition,  diffracted  rays  are  allowed 
that  originate  from  the  scattering  of  the  incident  wave  at  edges,  comers,  and  vertices. 
The  formulas  are  derived  on  the  basis  of  infinite  frequency  (A  — >  0)  which  implies  an 
electrically  large  target.  Among  the  problems  associated  with  GO  is  that  the  source  and 
the  observation  points  must  satisfy  Snell’s  Law  or  the  reflected  field  is  zero.  It  yields 
discontinuous  fields  at  shadow  boundaries  and  zero  fields  in  the  shadow  regions.  The 
Geometric  Theory  of  Diffraction  was  developed  to  take  into  account  the  edge-related 
effects  by  decomposing  the  total  field  at  the  observation  point  into  its  GO  and  diffracted 
components,  thereby  resulting  in  approximations  for  edge  diffracted  fields  [Refs.  1  and 
4]. 
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4.  Physical  Optics 

Physical  Optics  overcomes  the  problems  associated  with  the  infinities  of  ray 
tracing  for  flat  and  singly  curved  surfaces  by  approximating  the  induced  surface  current 
on  the  illuminated  portions  of  the  body  to  be  proportional  to  the  incident  magnetic  field 
intensity  and  integrating  them  to  obtain  the  scattered  far  field.  On  the  shadowed  portion 
of  the  target,  the  current  is  set  to  zero.  Thus,  the  current  is  approximated  as 

[in  x  Hi ,  for  the  illuminated  portion 
J  s  —  \  (A®) 

[0 ,  for  the  shadowed  portion 

where  Hi  is  the  incident  magnetic  field  intensity  at  the  surface  and  n  is  the  surface 
normal  unit  vector. 

The  approximation  in  equation  (2.8)  is  used  in  the  radiation  integrals  to  compute 
the  scattered  far  field  from  the  target.  PO  is  a  high-frequency  approximation  method  that 
gives  best  results  for  electrically  large  bodies  (L  >  1C )A) .  It  is  most  accurate  in  the 
specular  direction  and  gives  good  estimates  at  points  far  from  discontinuities  and  shadow 
boundaries.  However,  because  PO  abruptly  sets  the  current  to  zero  at  a  shadow 
boundary,  the  computed  field  values  at  wide  angles  and  in  the  shadow  regions  are 
inaccurate.  Also,  surface  waves,  multiple  reflections  and  edge  diffractions  are  not 
included  [Refs.  1,  2  and  4]. 
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III.  PHYSICAL  OPTICS  APPLIED  TO  TRIANGULAR  PATCHES 


A.  SCATTERING  FROM  FACETS  WITH  ARBITRARY  ORIENTATION 

The  geometry  of  a  single  facet  with  an  arbitrary  orientation  is  shown  in  Figure  2. 
The  facet  vertices  are  known  in  terms  of  their  cartesian  coordinates  ( xn ,yn,z„ )  for 

n  =  1...3.  It  is  assumed  that  the  vertices  are  labeled  in  a  right-hand  sense 
(counterclockwise)  so  that  the  thumb  of  the  right  hand  gives  the  outward  normal 
direction.  For  a  closed  body  that  surrounds  the  origin,  the  outward  normal  will  in  general 
point  away  from  the  origin.  The  position  vectors  to  the  nodes  are  given  by 
rn  =xnx  +  yny  +  znZ. 


Figure  2.  Arbitrarily  Oriented  Facet,  (a)  Arbitrarily  Oriented  Triangular  Facet, 

(b)  Facet  Sub-Areas  for  Natural  Coordinates 

The  edge  vectors  are  determined  from  the  node  coordinates  by 

\=r3-r2,  C3-1) 
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The  tip  of  edge  vector  n  originates  at  node  n.  The  outward  normal  is  obtained  by  taking 
the  cross  product  of  any  two  edge  vectors  in  a  right-handed  sense.  For  example, 

A  _ 

n-jg  =  nxx+nyy  +  nzz  (3.2) 

With  this  ordering  of  the  vertices,  the  outward  side  of  the  facet  is  denoted  as  the  front 
face  and  the  opposite  side  as  the  back  face.  If  a  plane  wave  is  incident  from  an  angle 
propagating  towards  the  origin,  then  its  propagation  vector  is  given  by 

A  t  * 

kt=~r  =  ~(x  sin  Ql  cos  fa  +  y  sin  6X  sin  fa  +  z  cos  0t )  =  -{xut  +  yv;-  +  zwt )  (3 .3) 

where  (m,  ,  v£- ,  wt )  are  the  x,  y  and  z  direction  cosines,  respectively,  and  r  the  radial  unit 
vector  from  the  origin  to  the  source  at  (0iyfa).  It  is  a  simple  matter  to  determine 
whether  the  front  face  is  illuminated  or  not  by  performing  the  test 

-ki»n>  0  (3.4) 

The  total  area  of  the  facet,  A,  can  be  computed  from  a  vector  relationship  between  edges 

A  =  |plx^|  (3.5) 

In  order  to  simplify  integration  over  the  surface  of  a  triangle,  normalized  area 
coordinates  (also  known  as  simplex  or  natural  coordinates)  are  used.  Consider  an 
integration  point  P  located  on  the  plate  at  (xp,yp,Zp) .  Its  position  vector  is 

rp  =xPx+yPy  +  zPz  (3.6) 

and 

P\=rp-r [,  p2=rP-  r2,  P2=rp-  r3  (3.7) 
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Note  that  for  the  monostatic  case 


g=h 

The  scattered  field  is  written  as 


Es  2  nr 


jk 


jkr 


V^e-Ej}\\e^dsf 


(3.15) 


(3.16) 


POLARIZATION  A 
FACTOR,  PF  ' - 


It  is  not  possible  to  obtain  an  exact  closed  form  result  for  the  integral.  However, 
given  that  the  incident  wavefront  is  plane  and  known  at  the  facet  vertices,  the  amplitude 
and  phase  at  interior  integration  points  can  be  found  by  interpolation.  Then  the  integrand 
can  be  expanded  using  a  Taylor  series,  and  each  term  integrated  to  give  a  closed  form 
result.  Usually  a  small  number  of  terms  (on  the  order  of  5)  will  give  a  sufficiently 
accurate  approximation  for  the  integral. 

Using  the  notation  of  [Ref.  5],  define 


h  =jjcc(TJ,Z)ejDc^)dr}dZ 


where 


Cc  (JJ,  €)  =  Cpr)  +  Cq£  +  CQ  = 
For  a  unit  amplitude  plane  wave 


Ei 


Similarly, 


Et  -1— »  Cp  -Cq  -0,  CQ  =1 


Dc(tj,£)  =  D  7] +  Dq£  +  D0 


(3.17) 


(3.18) 


(3.19) 


(3.20) 


where 


Dp  =  *[(*!  -x3)u  +  (yj  - y3 )v  +  (zi  - z3 )w] 
Dq  =  ^[(^2  -X3)u  +  (y2  -y3)v  +  (z2-  Z3  )w] 
Da  =  k[x3u  +  y3v  +  z3w] 
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(3.21) 


Now  the  integral  is  given  by 


\jdp 

Co 

_JDq 

Co 

C0 

1 

Dp{Dq-Dp)_ 

£ 

Dq(Dq-Dp)_ 

DqDp\ 

(3.22) 


To  avoid  the  numerical  problems  that  occur  near  the  singularities  (denominators  near 
zero),  the  following  special  cases  are  used: 

Case  1:1  D, 


h  = 


<Lt  and  |D?|  >Lt 
2  AejD°  ~ 

.  _  2Lt 


0 Dp) 

f -Co 

n\ 

[tt  +  1 

+  eJ  q  C0G(n,-Da ) 


Case  2: 


D, 


<  Lt  and  Dq 


<  L, 


le=2Ae^l 


C0(jDpY{jD„Y’ 


Case  3: 


D , 


n=0m=0  (m  +  n  +  2)\ 

>Lf  and  Dq  <Lt 

jD.JD ,  ”  ()DqT(  C„ 


lc=2Ae^e'  "  % 


n=0 


n\  [  n  + 1 


G(n  +  1,—Dp) 


Case  4:  \D. 


4  = 


—  Lf ,  Dq  >  Lt ,  and  Dq  Dp  <  Lt 


2 Ae»o  ~  O' Dp-jDqf 


jDq  n=  0  ra! 

The  function  G  is  defined  by 

1 

G(n,y)  =  j  sne^s  ds 


eiDqCn 


-C0G(n,Dq)  + - - 

^  n  +  1 


(3.23) 


(3.24) 


(3.25) 


(3.26) 


(3.27) 


0 

There  is  recursive  relationship  that  is  efficient  for  finding  G  for  higher  values  of  n 

n>\ 

jy 


(3.28) 
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and 


ejr  -1 
G(0,y)  =  - - t 

jy 


(3.29) 


C.  COORDINATE  TRANSFORMATIONS 

In  general,  the  local  facet  coordinate  system  will  not  be  aligned  with  the  global 
coordinate  system.  A  series  of  transformations  must  be  performed  in  order  to  obtain  the 
required  relationships  between  the  local  and  global  coordinate  variables  and  unit  vectors. 
Figure  3  shows  a  facet  with  arbitrary  orientation.  The  global  coordinates  are  (x,  y,  z)  and 
the  local  coordinates  are  (x",y",z") .  Figure  4  shows  that  the  z"  axis  is  aligned  with  the 
facet  normal.  Also,  the  figure  shows  a  right  triangle  facet  with  the  x"  and  y"  axes 
aligned  with  the  triangle  edges,  however,  the  transformation  can  be  applied  to  any 
arbitrary  facet  that  lies  in  the  x"  -  y"  plane. 


Figure  3.  Facet  with  Arbitrary  Orientation 
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Figure  4.  Facet  Rotation  Angles 


The  double-primed  variables  can  be  obtained  by  a  series  of  two  rotations  of  the 
global  coordinate  variables.  The  first  is  a  rotation  about  the  z  axis  by  an  angle  a  as 
shown  in  Figure  5.  The  resulting  coordinate  axes  are  referred  to  as  (x \y',z').  The 
rotation  angle  is  a-  atm2(ny,nx)  (atan2  is  the  MATLAB  arctangent  function  that  gives 

the  proper  quadrant).  The  second  rotation  is  by  an  angle  /?  =  cos_1(£  •  n)  about  the  y' 
axis,  illustrated  in  Figure  6. 


/ 


Figure  5.  First  Rotation  by  an  Angle  a  about  the  z  axis. 
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Figure  6.  Second  Rotation  by  an  Angle  (3  about  the  y'  axis. 


The  transformation  matrices  are  determined  from  geometry: 


r  /“ 

X 

cos  a 

sin  a 

O' 

X 

/ 

y 

= 

-sin  or 

cos  a 

0 

y 

/ 

_Z  _ 

0 

0 

1 

_z_ 

T' 


r 

X 

cos/? 

0 

-sin  / 3 

-  /' 
X 

y 

— 

0 

1 

0 

/ 

y 

rt 

Z  _ 

sin  ft 

0 

cos  P 

/ 

_Z  _ 

The  double  overbar  denotes  a  matrix.  Therefore, 


|V1 

X 

i - 

M 

*  * 
l _ 

_rr 

y 

_z_ 

The  same  transformation  matrix  holds  for  the  unit  vectors 


V 

X 

"  A  “ 

X 

*  rr 

y  * 

= rr 

y 

X . 

i_ 

(3.30) 


(3.31) 


(3.32) 


(3.33) 


20 


and  the  direction  cosines 


|V1 

u 

I 

*  V 

- 1 

=rr 

V 

w 

The  direction  cosines  are 


u  =  sin  #  cos  (p 
v  =  sin  (9  sin  ^ 
w  =  cos# 


(3.34) 


(3.35) 


and 

u  =  sm0  cos  (p 

v'1'  =  sin#'sin^',  (3.36) 

w/”  =  cos#'r 

These  quantities  can  be  either  incident  angles  (subscript  “z”)  or  observation  angles 
(subscript  “/?”)•  For  the  monostatic  case  both  are  equal  and  subscripts  are  not  required. 


D.  PHYSICAL  OPTICS  CURRENT 

The  components  of  the  surface  current  in  terms  of  the  local  cartesian  coordinate 
variables  are  from  [Ref.  6] 


and 


Eg.  cos  Q”  cos  f  Ej£  cos 
2 Rs  +  Za  cos  6 '  2 Rs  cos  Q"  +  Z0 


(3.37) 


_  ~ Eie.  cos#'1' sin  0*  Ej£  cos#" cos 0* 
**  ”  I  2 Rs  Tzo  cos#'  2 Rs  cos #'  +  Z0 


(3.38) 


where  Rs  is  the  surface  resistivity.  When  Rs  =  0  ( |T|  =  1) ,  the  surface  is  a  perfect 


electric  conductor.  As  Rs  — »  °°  the  surface  becomes  transparent  (|r|  — >  0) .  Defining  the 


transverse  magnetic  (TM)  and  transverse  electric  (TE)  reflection  coefficients  as 
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(3.39) 


rTM  =  r||  = 
rTE  =  ri  = 


~ZQ  cos  9" 
2RS  +  Z0  cos  6" 

ZlsL 

2  Rs  cos  6" +Z0 


allows  the  surface  current  components  to  be  expressed  as 


E,r  cosf^  ,  Eif  sinfr 

1  II  T  —  1  i 


/V/ 

COS# 


Jy  = 


Eie-  sin  <j>"  El0.  cos  f . 


■r!- 


cos^' 


The  spherical  polar  angles  are  known  from  the  direction  cosines 

sin^  =  V(«')2  +  (^)2 
cos^  =  Vl-(«')2+(v")2 


/  *  ' 

(  V 


tan  (/)"  =  arctan 

l 

In  the  bistatic  case,  the  incident  angles  are  used  in  the  above  formulas. 


(3.40) 


(3.41) 

(3-42) 


(3.43) 


E.  DIFFUSE  CALCULATION 

Models  may  have  a  rough  surface  due  to  irregularities  or  material  imperfection. 
Surface  roughness  introduces  path  differences  relative  to  a  smooth  surface  and  thus 
contribute  a  diffuse  scattering  component  to  the  RCS  that  is  added  to  the  coherent 
component.  The  random  deviation  of  the  surface  results  in  diffuse  scattering  which 
yields  a  relatively  angle  dependent  term  to  the  RCS.  First  we  assume  that  the  surface 
deviation  resemble  a  Gaussian  shape  as  shown  in  Figure  7.  The  extent  of  the  Gaussian 
shape  is  characterized  by  a  correlation  distance  C  (in meters),  which  is  the  average 
distance  at  which  deviations  become  uncorrelated.  A  large  correlation  distance  implies  a 


22 


slowly  varying  surface  error,  whereas  a  small  correlation  distance  is  associated  with  a 

2 

rapidly  varying  error.  The  variance  of  the  irregularities  is  denoted  as  8  (or  standard 
deviation  as  8)  as  illustrated  in  Figure  8. 


Figure  7.  Gaussian  Approximation  of  Rough  Surface  From  [Ref.  1] 


— *■  * 


Figure  8.  Correlation  Distance,  (a)  Large  Correlation  Distance, 
(b)  Small  Correlation  Distance  From  [Ref.  1] 


In  the  PO  calculation,  surface  roughness  is  approximated  by  defining  values  for 
correlation  distance  and  standard  deviation.  The  average  value  of  the  RCS  is 
approximately  [Ref.  7] 


a 


ave 


4  ttA1 
~ 


cos2  6  e 


-4 k'-s- 


4k282nC2  f 
— expj- 


(3.44) 
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where  A  is  the  area  of  the  facet  and  6  is  measured  from  the  normal.  This  value 
simply  added  to  the  specular  term 


F.  THE  SCATTERED  FIELD  OF  A  COLLECTION  OF  FACETS 

An  arbitrary  target  can  be  approximated  by  a  collection  of  triangular  facets. 
Equation  3.16  must  be  evaluated  for  all  facets,  and  the  polarization  factors  and  coordinate 
systems  may  differ  for  each  facet.  However,  the  integral  over  all  facets  will  have  the 
form  of  equation  3.17.  To  obtain  the  scattered  field  of  the  array  of  facets  the  following 
procedure  is  used: 

1.  The  incident  field  is  specified  in  the  global  coordinates:  angle  of  incidence 
and  polarization  (Ei0,Ei(p). 

2.  For  each  facet  (indexed  from  c  =  l,2,...,Nc)  the  incident  field  is 

transformed  to  the  local  facet  coordinate  system.  Also  the  angle  of  incidence  is 
determined  in  the  local  facet  coordinate  system. 

3.  The  physical  optics  current  is  computed  using  the  incident  magnetic  field,  and 
the  radiation  integral  evaluated  over  the  facet  surface  to  obtain  the  scattered  field 
components  (Esq*,ES0')  . 

4.  The  scattered  field  components  are  transformed  back  to  the  global  coordinate 
system,  thereby  obtaining  (EsQ,Es(p)  for  each  facet. 

5.  All  of  the  Nc  transformed  scattered  field  components  are  superimposed 

(vector  summed)  to  obtain  the  total  scattered  field  in  the  observation  direction. 

6.  The  RCS  is  computed  from  the  total  scattered  field  from  equation  (2.3). 

A  flow  diagram  of  the  procedure  for  the  monostatic  calculation  is  shown  in  Figure  9.  For 
a  bistatic  calculation,  the  current  does  not  have  to  be  recomputed  for  each  observation 
angle  because  the  incident  angle,  and  hence  current,  is  fixed.  [Ref.  8] 
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Figure  9.  Flow  Chart  for  Computing  Monostatic  RCS  of  an  Array  of  Facets 


IV.  THE  POFACETS  GRAPHICAL  USER  INTERFACE 


POFACETS  is  a  graphical  user  interface  (GUI)  radar  cross  section  (RCS) 
prediction  code  based  on  the  physical  optics  (PO)  approximation  and  implemented  in 
MATLAB.  It  is  accurate  when  the  target  is  large  in  terms  of  wavelength.  Scattering 
objects  are  approximated  by  arrays  of  triangles  (facets)  and  superposition  is  used  to 
compute  the  total  RCS  of  the  model  object  as  described  in  Chapter  III.  Multiple 
reflections  are  not  included.  The  GUI  consists  of  three  modules:  (1)  Design  Model ,  (2) 
Calculate  Monostatic  RCS,  and  (3)  Calculate  Bistatic  RCS. 

Design  Model  provides  a  convenient  way  of  modeling  a  complex  object  or  loading 
and  editing  a  previously  designed  model.  The  geometry  of  the  model  is  defined  using 
two  sets  of  data.  The  first  set  defines  the  x,y,z  coordinates  of  the  vertices.  The  second 
set  of  data  defines  the  component  triangular  facets  of  the  object,  their  nodes,  and 
illumination  and  resistivity  characteristics. 

Calculate  Monostatic  RCS  computes  the  total  radar  cross  section  of  the  target  using 
a  monostatic  radar  (co-located  transmitter  and  receiver).  The  scattered  field  of  each 
triangle  is  computed  as  if  it  was  isolated  and  other  triangles  were  not  present.  Multiple 
reflections,  edge  diffraction  and  surface  waves  are  not  included.  Shadowing  is 
approximately  included  by  considering  a  facet  to  be  completely  illuminated  or 
completely  shadowed  by  the  incident  wave.  Options  for  defining  the  roughness  of  the 
surface  are  available. 
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Calculate  Bistatic  RCS  computes  the  total  radar  cross  section  of  the  target  using  a 
bistatic  radar  (transmitter  and  receiver  at  different  locations).  The  incidence  direction  is 
fixed  and  the  observation  direction  changes. 

To  run  the  POFACETS  GUI,  the  user  must  first  run  MATLAB  and  set  the  path  or 
directory  to  where  the  POFACETS  code  is  stored.  This  may  be  done  in  two  ways.  The 
first  is  by  using  the  File  Menu  and  selecting  Set  Path.  This  displays  the  Path  Browser 
window  and  from  here  click  on  the  Browse  button.  In  the  Browse  for  Folder  window, 
select  the  folder  where  POFACETS  is  located  then  click  the  OK  button.  The  second 
method  is  by  using  the  change  directory  command  at  the  MATLAB  command  prompt  to 
where  the  code  is  stored  such  as  cd.  c:\pofacets.  After  setting  the  path  or  directory,  typing 
pofacets  at  the  command  line  will  execute  the  GUI  as  shown  in  Figure  10.  A  brief 
explanation  of  the  POFACETS  GUI  may  be  obtained  if  the  user  clicks  on  the  Help 
button,  which  displays  the  POFACETS  Help  GUI  as  shown  in  Figure  1 1 . 
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-J>  RCS  PREDICTION  USING  PHYSICAL  OPTICS  APPROXIMATION 


EkaoE.  Garrido  Jr.  and  David  C.  Jean 


Create  a  new  target  geometry  or  load; 
and  edit  a  previously  designed  model. 


Design  Model 


Compute  fee  monostatic  radar  crass 
section  of  fee  modeLl 


Calculate  Monostatic  RCS 


Figure  10.  POFACETS  Main  GUI 
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PO FACETS  Help 


I POFACETS  is  a  MATLAB  graphical  user  interface  (GO I)  radar  cross 


section  prediction  (JR.CS)  code  based  on  the  physical  optics  (PO) 
approximation  method.  It  is  accurate  for  relatively  high  frequencies. 
Scattering  objects  are  approximated  by  arrays  of  triangles  (facets) 
and  superposition  is  used  to  compute  the  total  RCS  of  the  model 
object 

Design  Model  provides  a  convenient  way  of  modeling  a  complex  object 
or  to  load  and/or  edit  a  previously  designed  and  saved  model  The 
geometry  of  the  model  is  defined  using  two  sets  of  data.  The  first 
set  defines  the  xyz  coordinates  of  the  vertices.  The  second  set  of 
data  defines  the  component  triangular  meets  of  the  object,  its  nodes, 
illumination  and  resistivity  characteristics. 

Calculate  Monostatic  1RC$  computes  the  total  radar  cross  section  of 
the  target  using  a  monostatic  radar  (co-located  transmitter  and  receiver). 
The  scattered  field  of  each  triangle  is  computed  as  if  it  were  isolated 
and  other  triangles  were  not  present.  Multiple  reflections,  edge 
diffraction  and  surface  waves,  are  not  included  Shadowing  is  approximate! 


Figure  11.  POFACETS  Help  GUI 


A*  DESIGNING  AND  PLOTTING  A  MODEL 

To  illustrate  the  Design  Model  module  of  the  GUI,  we  will  design  the  AIM-9 
Sidewinder  missile  shown  below  in  Figure  12  by  modeling  it  with  the  set  of  coordinates 
in  Table  1  and  facet  configuration  in  Table  2. 


Figure  12.  AIM-9  Sidewinder  from  [Ref.  9] 
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Vertex  # 

X  (m) 

Y(m) 

Z(m) 

Y(m) 

Z(m) 

1 

0.000 

0.000 

0.00 

20 

-0.250 

0.065 

0.150 

2 

0.130 

0.000 

0.00 

21 

-0.250 

0.065 

0.000 

0.130 

0.130 

0.00 

22 

0.065 

0.000 

0.000 

0.000 

0.130 

0.00 

23 

0.065 

0.000 

0.200 

0.000 

0.000 

2.57 

24 

0.065 

-0.250 

0.150 

6 

0.130 

0.000 

2.57 

25 

0.065 

-0.250 

0.000 

0.130 

0.130 

2.57 

26" 

0.130 

0.065 

2.200 

8 

0.000 

0.130 

2.57 

27 

0.130 

0.065 

2.400 

0.065 

0.065 

2.87 

28 

0.380 

0.065 

2.200 

10 

0.130 

0.065 

0.00 

29 

2.200 

11 

0.130 

0.065 

020 

30 

0.130 

0.065 

2.400 

12 

0.380 

0.065 

015 

31 

0.065 

0.380 

2.200 

13 

0.380 

0.065 

0.00 

32 

0.000 

0.065 

2.200 

14 

0.065 

2.400 

15 

0.065 

0.130 

0.200 

34 

-0.380 

0.065 

2.200 

16 

0.065 

0.385 

0.150 

35 

0.065 

0.000 

2.200 

17 

0.065 

0.385 

0.000 

36 

0.065 

0.000 

2.400 

18 

0.000 

0.065 

0.000 

37 

0.065 

-0.380 

2.200 

19 

0.000 

0.065 

0.200 

Table  1.  AIM-9  Sidewinder  Coordinate  Set 
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Facet  # 

Node  1 

Node  2 

Node  3 

Illumination 

Normalized 

Resistivity 

1 

3 

2 

1 

1 

1 

2 

1 

4 

3 

i 

i 

3 

3 

7 

6 

i 

i 

4 

6 

2 

3 

i 

i 

5 

4 

8 

7 

1 

1 

6 

7 

3 

4 

i 

i 

7 

4 

1 

5 

i 

i 

8 

5 

8 

4 

i 

i 

9 

2 

6 

5 

1 

1 

10 

5 

1 

2 

1 

1 

11 

2 

9 

6~ 

1 

1 

12 

8 

9 

7 

i 

1 

13 

5 

9 

8 

1 

1 

14 

5 

6 

9 

1 

1 

15 

10 

11 

12 

0 

i 

16 

12 

13 

10 

0 

1 

17 

16 

15 

14 

0 

i 

18 

14 

17 

16 

0 

1 

19 

19 

18 

6 

i 

20 

18 

21 

20 

0 

1 

21 

22 

23 

24 

0 

1 

22 

24 

25 

22 

0 

1 

23 

26 

27 

28 

0 

i 

24 

29 

31 

30 

0 

i 

25 

32 

34 

33 

0 

i 

26 

35 

36 

37 

0 

1 

Table  2.  AIM-9  Sidewinder  Facet  Configuration 
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To  design  a  new  object,  the  user  clicks  on  the  Design  Model  button  of  Figure  10. 
This  displays  the  Design  Model  Coordinates  GUI  shown  in  Figure  13. 


Figure  13.  Design  Model  Coordinates  GUI 


The  Number  of  Vertices  is  initially  disabled  but  selecting  File  j  New  in  the  menu 
will  enable  this  box  to  accept  an  entry.  Since  the  Sidewinder  example  has  37  vertices,  we 
enter  this  number  in  the  appropriate  box  and  the  GUI  generates  the  corresponding 
number  of  editable  text  boxes  (all  initialized  to  zero)  to  accommodate  the  Cartesian 
coordinates  of  the  model  as  well  as  enable  the  Apply  button  to  check  all  the  individual 
coordinates  later.  Entering  the  coordinates  correctly  and  clicking  on  the  Apply  button 
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displays  the  Coordinate  Status  GUI  (Figure  14)  and  after  clicking  the  OK  button,  enables 
the  Design  Facets  button  and  the  active  GUI  should  now  look  like  Figure  15. 


Coordinate  Status 


mm 

Coordinates  valid,  press  Design  Facets  button . 
OK  | 


Figure  14.  Coordinate  Status  GUI 


DESIGN  MODEL  Coordinates 


Fife  .'E-EiE' 

Number  of  Vertices:  J  37  1 .j— 7 — 

'  4  VERTEX  COORDINATES  (m>  SCaRK  1  1 

vF  . . X  .  Y  2  v#  X  Y  :  W 

1 . I.  1 . . 0  M  . o  1  ,  .  0} 

1. 21  jv| . _ .  -0.25  j  . . .  0.0651  C  . 0 

1  2  1  0.13  1  0  I  0 

Eli -E. ' 0.035  C  0  f~  0 

!  3  ,1  0.13  j  0.13  j  0 

I  23  lE .  0165  j . .  0/4.....  . 0.2  .S 

1  i/  I .  0:|  ....  0.13  1  _ 0  * 

. . . 

1,1.1...  0  1  0  1  257 

1:25  m  0,065  1  *0.25  1  0 

■  El.  1.,  . 0.13  j  . 0:4  257 

\.m  1 _ 111;  L _ 0.065  l . . 22 

1  ?  1  0.13  J  013  J  257 

?  fi  i  ,  0  j  0.13'  j  257 

rf*  i .. .  q-131;  i .. ."  0.06=  j:..::e4. 

Jlp  1 . “ . '  038  f~  ,  0,065  r  .  2.2 

1  S.  | . ....0,035  1 _ 0065  j _ 287  ; 

.  Oil  1  0.13  j  0065  1  0 

JW.  1  0.065 ;  j~  ....0.13  t  2.2 

Tm  1 0055,.  j  ...  ...0.13  .  J  24 

JW  f  0.055  1  0.38  1  22 

j  -2  1  0.38  J  0.065  :  •  1  0.15 

132  ]..... . . 0  J  . 0065  J _ _ 22 

1  33  1... . . . 0j4, . ..1165  r  . 2.40. 

I  G.0S5  j  0.13  1  0 

1  34  )  -038  ■  j  0065  1  22 

El  1 . . .0.065  1 _ .0,13.  L _ 1,2 

111:  j . 0,065  i  .  0  i  22 

1  11;  1 . .  0065.  j . 0.385.  j . 01.5 

1.23.4 . . 0,065  J .  0  j .  24 

1  17!  |  0,065  .,  El  .  0.385  1  0 

::  OE  1  o  j  0.065  jE  0 

JIT  j  0065  r  -0.38  r  2.2 

1  19]  j  0  1  0.065  1  0.2 

[20  J  -025,.  I  0065  |  0.15 

;.E  Design  Facets  Apply  : 

Close  Help 

Figure  15.  Design  Model  Coordinates  GUI  with  Vertices  of  the  Sidewinder 
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*  DESIGN  MODEL  Facets 


is  enabled.  The  nodes  and  illumination  values  are  then  encoded  using  the  data  in  Table  2. 
The  illumination  parameter  specifies  whether  both  sides  of  a  facet  are  illuminated  by  the 
incident  wave.  If  Hum  =  0  both  sides  of  the  plate  are  illuminated  (e.g.,  a  fin  or  a  wing); 
if  Hum  =  1  only  the  outer  face  is  illuminated  (e.g.,  a  fuselage).  The  normalized  resistivity 
Rs ! 377  is  used  to  define  the  resistivity  characteristics  of  the  model  surface.  For 
example,  to  conveniently  define  the  entire  missile  skin  to  have  an  impedance  of  377 
ohms,  uniform  resistivity  is  set  to  1  and  it  will  be  noted  that  the  resistivity  values  for  all 
the  facets  are  accordingly  set  in  the  Rs  column.  Pressing  the  Apply  button  displays  the 
Facet  Configuration  Status  GUI  (Figure  17)  if  there  are  no  errors  and  subsequently 
clicking  the  OK  button  enables  the  View  Model  button  and  the  active  GUI  should  now 
look  like  Figure  18. 


**#  Facet  Configuration  Status 


Ena 


All  nodes,  illumination,  end  resistivity  values  valid.  Click  View  Model 
Sutton.  .■  '  •: 


OK 


Figure  17.  Facet  Configuration  Status  GUI 


To  visualize  the  model  comprised  of  its  triangular  components,  the  plot  may  be 
seen  by  clicking  the  View  Model  button,  which  displays  the  Model  Geometry  GUI  as 
shown  in  Figure  19. 
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Figure  18.  Design  Model  Facets  GUI  with  Facet  Configuration  of  the  Sidewinder 


Figure  19.  Model  Geometry  GUI 
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1.  Help  GUIs 

A  number  of  Help  GUIs  are  available  to  provide  the  user  valuable  assistance  in 
working  with  the  Design  Model  module.  The  three  Help  GUIs  may  be  viewed  by 
clicking  on  the  Help  buttons  of  the  Design  Model  Coordinates,  Design  Model  Facets,  and 
the  Model  Geometry  GUIs  to  display  Figures  20,  21  and  22,  respectively. 


Figure  20.  Design  Model  Coordinates  Help  GUI 
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-#  Design  Model  Facets  Help 

j  Design  Model  Facets  conjures  the  geometry  of  the  object  by  defining 


an  array  of  triangular  components  or  facets  comprising  the  model  The 
facets  are  described  using  a  facet  number,  which  by  default  is 
sequential.  Each  facet  is  defined  by  three  nodes  or  vertices,  and  its 
illumination  and  resistivity  characteristics. 

Nodes  -  Entries  in  the  nodes  (no del,  node 2,  and  node 3)  represent  the 
vertex  numbers  defining  the  facet  The  nodes  comprising  the  facet 
should  be  ordered  n  a  right-hand  or  counter  clockwise  sense  where 
the  normal  to  the  facet  is  the  direction  of  illumination. 

Bumination  (Sum)  -  The  ffiuminarion  value  is  0  if  the  facet  is 
allowed  to  be  illuminated  from  both  sides  and  1  if  the  meet  will 
only  be  illuminated  from  one  side. 


Resistivity  (Rs)  -  The  resistivity  of  the  meet  is  normalized  to  the 
free  space  value  of  37?  ohms.  If  surface  Rs  =  0,  then  the  facet  is  a 
perfect  electric  conductor  (PEC);  when  Rs  is  infinity  (in practice, 
several  thousand)  the  facet  is  transparent 


Close 


Figure  21.  Design  Model  Facets  Help  GUI 


*  Model  Geometry  Help 


The  plot  of  the  Model  Geometry  is  a  scaled  version  of  the  design 


defined  by  the  coordinates  and  the  facets  configuration. 

Label  Vertices  -  is  a  function  that  plots  the  model  geometry  showing 
the  vertex  numbers  as  defined  in  the  Design  Target  Model-Coordinates 
GDI 

Labe!  Faces  -  is  a  function  that  plots  the  model  geometry  showing  the 
facet  or  triangle  number  as  defined  by  the  nodes  in  the  Design  Target 
Model  -  Facets. 

Grid  On/Off-  is  a  function  that  plots  or  clears  the  axis  area  of  the 
model  geometry  figure  of  grid  fines  to  enhance  visualization  of  the 
model. 

Print  Files  -  this  function  prints  a  hard  copy  of  the  c  o  ordinate  s/£ac  et  s 
data  and/or  the  model  geometry  as  selected 

Close  -  doses  the  GDI  figure. 


l'}A' 


vVy; 


wv. 


.Close 


Figure  22.  Model  Geometry  Help  GUI 
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2. 


Status  GUIs 


Displaying  several  Status  GUIs  have  been  incorporated  in  the  code  to  draw  user 
attention  to  specific  part  of  the  GUI  requiring  action.  For  the  Design  Model  Coordinates 
GUI,  an  error  in  the  Number  of  Vertices  entered  will  display  the  status  shown  in  Figure 
--*■  Likewise,  since  the  Design  Model  Coordinates  GUI  can  initially  accommodate  a 
maximum  of  40  vertices  only,  exceeding  this  limit  will  display  Figure  24.  Meanwhile, 
encoding  an  invalid  Scale  value  will  show'  Figure  25. 


I  Number  of  Vertices  Status 

sssi 

Please  check  the  numb 

.  OK  j 

er  of  vertices  entered!  u 

Figure  23.  Number  of  Vertices  Status  GUI 


Maximum  Number  of  Verices 


as m 


A  Design  can  handle  up  to  40  vertices  only! 


:  OK 


Figure  24.  Maximum  Number  of  Vertices  Status  GUI 


I  Scale  Status 

mmmw 

O 

:‘lease  check  scale  of  vertices  entered.  i 

;  OK  | 

Figure  25.  Scale  Status  GUI 
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The  user  is  advised  of  erroneous  entries  on  the  coordinates  with  the  display  of 
Figures  26a  to  c,  while  a  valid  set  of  coordinates  is  ascertained  by  Figure  26d. 


(d) 


Figure  26.  Coordinate  Status  GUIs,  (a)  Invalid  Coordinate,  (b)  Uncorrected  Erroneous 
Coordinate,  (c)  Duplicate  Set  of  Coordinates,  (d)  Valid  Set  of  Coordinates 


For  the  Design  Model  Facets  GUI,  an  error  in  the  Number  of  Facets  entered  will 
display  Figure  27.  Likewise,  since  the  GUI  can  initially  accommodate  a  maximum  of  40 
facets  only,  exceeding  this  limit  will  display  Figure  28.  Meanwhile,  encoding  an  invalid 
resistivity  value  for  a  particular  facet  will  show  Figure  29. 


Figure  29.  Uniform  Resistivity  Status  GUI 


The  user  is  advised  of  erroneous  entries  in  the  configuration  of  the  facet(s)  with 
messages  shown  in  Figures  30a  to  f  and  encoding  valid  parameters  is  ascertained  by 
Figure  30g. 


42 


Figure  30. 


(a) 


Node  Status 


0 


I  ■■■  3Ifl 

Duplicate  nodes  for  the  same  facet  number! 
OK  i 


(b) 


Node  Status 


□m 


Error  in  node(s),  illumination  or  resistivity  entered!  Correct  entries 
where  9993  is  displayed. 


OK 


(C) 


(d) 


Facet  Configuration  Status  GUIs,  (a)  Invalid  Node,  (b)  Duplicate  Nodes, 
(c)  Unconnected  Erroneous  Node,  (d)  Duplicate  Set  of  Nodes 
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(§) 

Figure  30  (continued).  Facet  Configuration  Status  GUIs,  (e)  Invalid  Illumination  Value, 
(f)  Invalid  Facet  Resistivity,  (g)  Valid  Facet  Configuration 


3.  Model  Geometry  GUI  Functionalities 

The  Model  Geometry  GUI  shown  in  Figure  19  provides  some  functionalities  to 
enhance  user  visualization  of  the  model  object  and  to  verify  correctness  of  node  and  facet 
designation.  These  include  plotting  the  vertex/facet  numbers  with  the  object  or  the  grid 
lines  by  selecting  any  or  all  of  the  following  checkboxes:  Label  Vertices,  Label  Faces, 
and  Grid,  respectively.  Selecting  the  Coordinates/Facets  and  Model  Geometry 
checkboxes  and  clicking  on  the  Print  button  will  generate  a  hard  copy  of  the  coordinates 
and  facets  data. 
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4.  Saving  the  Model  Files 

The  model’s  file  is  composed  of  two  sets  of  data,  coordinate  set  and  facet 
configuration.  However,  to  save  a  model  only  the  filename  is  required.  To  save  the 
Sidewinder  example,  click  on  File  |  Save  As ...  from  the  menu  and  this  displays  the  New 
File  GUI  shown  in  Figure  31.  We  enter  the  filename  Sidewinder  and  click  on  the  Save 
button.  The  program  saves  the  files  and  confirms  the  action  to  the  user  by  displaying  the 
File  Saved  GUI  shown  in  Figure  32.  Inadvertently  specifying  a  filename  that  already 
exists  in  the  list  of  models  will  prompt  the  user  to  overwrite  said  file  with  the  Replace 
File  GUI  shown  in  Figure  33. 


£ 

.. 

is?,"": 

BOX 

1 

WAIL 

— I 

PLATE 

BRICK 

CHUCK _  .  _ 

zl 

File  name : 

Save  | 

Sidewinder] 

Cancel 

Figure  31.  New  File  GUT 


Figure  32.  File  Saved  GUI 
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c:^o1aeels^odels^LONG_BOXaIready  exists. 
Do  you  want  to  replace  It?  \ 


Figure  33.  Replace  File  GUI 


5.  Editing  an  Existing  Model 

To  edit  an  existing  model,  its  set  of  data  needs  to  be  retrieved  and  displayed 
before  editing  or  scaling  can  proceed.  This  is  accomplished  by  clicking  on  the  File  | 
Open  menu  of  the  Design  Model  Coordinates  GUI  (Figure  10)  which  displays  the  Open 
Model  File  GUI  shown  in  Figure  34.  Suppose  we  would  like  to  edit  the  VTAIL  model. 
We  select  it  among  the  listed  file  names  and  click  the  Open  button.  This  action  performs 
a  series  of  operations  that  will  display  its  data  set  and  the  model  plot  as  shown  in  Figure 
35. 


Figure  34.  Open  Model  File  GUI 
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Figure  35.  Displayed  VI AIL  File 


At  this  point,  the  user  may  click  on  any  of  the  three  overlapping  GUIs  to  make  it 
active  and  edit  the  pertinent  data.  The  different  Status  GUIs  mentioned  earlier  will  guide 
the  user  throughout  the  use  of  the  interface. 

B.  CALCULATING  THE  MONOSTATIC  RCS 
1.  Smooth  PEC  Model 

To  compute  the  Monostatic  RCS  of  a  model  object,  the  user  clicks  on  the 
Calculate  Monostatic  RCS  button  of  the  POFACETS  Main  GUI  (Figure  10).  This  clears 
the  screen  of  other  GUIs  (except  the  Main  GUI)  and  brings  up  the  Calculate  Monostatic 
RCS  GUI  with  default  values  and  parameters  shown  in  Figure  36. 
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Figure  36.  Calculate  Monostatic  RCS  GUI 


For  this  case,  we  will  calculate  the  monostatic  RCS  of  a  closed  box  model  with  a 
smooth  surface,  measuring  100  meters  on  each  side  and  which  is  made  of  perfect  electric 
conductor  (PEC)  material.  Since  this  is  a  closed  form  object,  only  the  outer  face  will  be 
illuminated  and  all  the  facets  will  have  a  uniform  illumination  of  1.  We  further  assume 
that  the  incident  wave  is  0 -polarized  (TM-  z ),  the  frequency  of  interest  15  MHz  and  we 
would  like  the  RCS  measured  at  o  —  0  plane  while  varying  6  from  0  to  360  degrees  to 
get  a  view  of  the  model  RCS  around  its  sides. 

Setting  up  the  GUI  with  the  parameters  mentioned  should  result  in  the  Calculate 
Monostatic  RCS  GUI  that  looks  like  Figure  37  and  the  RCS  plot  shown  in  Figure  38. 
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Two  computational  parameters  that  may  be  varied  are  the  Taylor  Series  Region  and  the 
Number  of  Terms  in  the  Taylor  Series.  See  [Ref.  5]  for  details  on  how  these  parameters 
affect  the  accuracy  and  computation  time.  In  general,  the  default  values  give  good 
results  with  relatively  fast  run  times.  Also,  a  two-dimensional  contour  can  be  generated 
by  setting  6  values  from  0  to  90,  the  <b  values  from  0  to  360,  increment  angle  to  2 
degrees,  and  the  resulting  computation  of  the  monostatic  RCS  is  shown  in  Figure  39. 


¥1  CALCULATE  MONOSTATiC  RCS 


Model  Name: 


#i 

BRICK 

CHUCK 

DART 

; :::: 

LONG  BOX 

PLATE 

Sidewinder 

UTA2 

w 

VTAIL 

. . .3 

Lead  File 


Surface  Roughness 

Coirdaiion  Dist.(m).j  0 
Standard  Dew.  (m):  I  0 


Theta 

Phi 

Starting  Angie  >■  j  0 

deg 

1  0 

deg 

Ending  .Angle  |  360 

"deg  T 

1  0 

deg 

increment  Angle  2 

;  de  g ' 

i 

deg 

5  Computational  Parameters 

Taylor  Series 

Length  of  Region:  J  C.Cd 
Number cfTerms:  "S  5 


Incident  Polarization 
j  Theta  (TM-z)  zl 


Frequency  j  0.015 

.■.*.4;:;.:. .  GHz  ' 


Calculate  RCS 


Print 


Close 


Figure  37.  Calculate  Monostatic  RCS  GUI  for  Smooth  PEC  Box  Model 
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Figure  38.  Monostatic  RCS  of  Smooth  PEC  Box  Model 


*  MONOSTATIC  RCS 


Figure  39.  Contour  Plot  of  Smooth  PEC  Box  Model 
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RCS-theta  is  specified  as  Gee  (#  polarized  incident  wave  and  #  polarized  receiver). 
This  is  the  co-polarized  component.  RCS-phi  is  o^e  (0  polarized  incident  wave  and  <p 
polarized  receiver).  The  latter  is  the  cross-polarized  component.  In  this  particular 
example,  it  happens  to  be  at  a  level  below  the  specific  threshold  in  the  code. 

2.  Model  with  Surface  Impedance 

An  object  coated  with  a  radar  absorbing  material  or  having  a  skin  made  of 
resistive  material  may  be  approximated  by  defining  its  surface  to  have  an  impedance 
value.  Equations  3.37  and  3.38  describes  the  components  of  the  surface  current  at  the 
target  which  is  inversely  proportional  to  the  resistivity  of  the  surface.  Using  the  box 
model  as  an  example  and  defining  its  surface  to  have  an  impedance  three  times  the  value 
of  free  space  (i.e.  uniform  resistivity  is  set  to  3),  we  see  that  the  monostatic  RCS 
calculation  (Figure  40)  is  15dB  down  in  the  specular  return  and  from  15  to  20dB  at  the 
other  various  angles  compared  to  the  PEC  model  shown  in  Figure  38. 

3.  Rough  Surface  Model 

Meanwhile,  to  observe  the  effect  of  the  diffused  scattering  component  of  a  rough 
surface,  the  box  model  is  defined  to  have  a  correlation  distance  of  100  meters  (or 
equivalently  5 A.)  and  a  standard  deviation  about  the  mean  value  of  2  meters  (that  is, 
0. 1A  for  a  frequency  of  15  MHz ).  Without  varying  the  other  parameters,  the  monostatic 
RCS  plot  in  Figure  41  shows  a  decline  of  about  8  dB  in  the  specular  direction  and  an 
increase  in  the  RCS  at  all  the  other  angles  due  to  the  diffuse  scattering  effect  compared  to 
the  PEC  model  results  in  Figure  38. 
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Figure  40.  Monostatic  RCS  of  Smooth  Box  Model  with  Surface  Impedance 


Figure  41.  Monostatic  RCS  of  Box  Model  with  Rough  Surface 
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4. 


Calculate  Monostatic  RCS  GUIs 


Several  GUIs  were  designed  to  enable  the  user  to  pinpoint  specific  parameters 
requiring  attention.  Errors  in  encoding  values  for  the  theta  and  phi  angles  display  Figures 
42a  to  c,  while  setting  the  same  starting  and  ending  angles  advise  the  user  that  a  theta  or 
phi-cut  is  set  as  shown  in  Figure  43. 


(c) 

Figure  42.  Angle  Status  GUIs,  (a)  Invalid  Theta  Starting  Angle,  (b)  Invalid  Theta  Ending 

Angle,  (c)  Invalid  Theta  Increment  Angle 
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Figure  43.  Theta-Cui  Set  GUI 


Other  GUIs  include  notifying  the  user  of  the  need  to  correct  entries  for  the 
Frequency,  Wavelength,  and  the  Length  of  the  Taylor  Series  Region  shown  in  Figures  44, 
45,  and  46,  respectively. 


Figure  45.  Wavelength  Status  GUI 


Figure  46.  Taylor  Series  Region  GUI 
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In  addition,  some  GUIs  have  been  pre-defined  to  warn  the  user  if  the  parameter 
entered  will  take  some  time  to  generate  the  results  as  in  Figure  47  or  suggest  a  range  of 
values  such  as  Figures  48  and  49. 


Taylor  Series  Terms 


:.lslx| 


Monosiatic  RCS  computation  for  number  of  terms  greater  then  10  may  take 
some tsime. : ' ■? : A-  V' 


OK 


Figure  47.  Taylor  Series  Terms  GUI 


U]  Correlation  Distance  Status 


T 

* 


Suggest  that  this  value  be  smaller  than  the  wavelength. 

OK 


Figure  48.  Correlation  Distance  Status  GUI 


Standard  Deviation  Status 


nsm 


A\  Value  entered  Is  valid  bii  extremely  high. 


K  i:  OK 


Figure  49.  Standard  Deviation  Status  GUI 


It  will  be  evident  to  the  user  that  the  Calculate  RCS  button  is  initially  disabled. 

Selecting  a  model  from  the  list  box  and  clicking  the  Load  File  button  will  confirm  to  the 

user  the  model  selected  and  give  a  reminder  to  set  all  the  parameters  as  shown  in  Figure 

50.  Clicking  the  OK  button  will  enable  the  Calculate  RCS  button.  Finally,  clicking  the 

Calculate  RCS  button  informs  the  user  that  the  GUI  is  calculating  the  monostatic  RCS  of 
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the  model  selected  as  shown  in  Figure  51  below.  The  Calculate  Monostatic  RCS  Help 
GUI  showm  in  Figure  52  can  also  provide  some  valuable  information. 


Figure  51.  Computing  Monostatic  RCS  GUT 


Calculate  Monostatic  RCS  Help 


f  Calculate  Monostatic  RCS  computes  the  total  radar  cross  section  of  a 


■  model  using  a  monostatic  radar  (collocated  transmitter  and  receiver). 
Scattering  objects  are  approximated  by  arrays  of  triangles  (facets)  and 
trie  scattered  field  of  each  triangle  is  confuted  as  if  it  were  isolated 
and  other  triangles  were  not  present.  Multiple  reflections,  edge 
di fraction  and  surface  waves  are  not  included.  Shadowing  is 
approximately  included  by  considering  a  facet  to  be  completely 
:  illuminated  or  completely  shadowed  by  trie  incident  wave. 


Start  and  Ending  Angles  -  The  illuminated  surface  of  a  model  is  described 
by  defining  the  starting  and  ending  angles  (theta  and  phi  values)  in  degrees 
:  of  the  model  with  respect  to  trie  monostatic  radar.  Angles  are  measured  in 
degrees  and  entering  a  phi  start  and  phi  stop  with  equal  values  results  in  a 
phi-cut.  Similarly,  entering  a  theta  start  and  theta  stop  with  equal  values 
:  results  in  a  theta-cut. 

:  Increment  -  defines  the  number  of  degrees  (or  fraction  of  a  degree)  increase 
in  step  for  the  start  and  ending  angles  to  be  used  in  calculating  the  mono  static 
radar  cross  section  of  the  model. _ _ 


Figure  52.  Calculate  Monostatic  RCS  Help  GUI 
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C.  CALCULATING  THE  BISTATIC  RCS 

The  Calculate  Bistatic  RCS  GUI  (Figure  53)  is  very  similar  to  the  Calculate 
Monostatic  RCS  GUI  (Figure  36)  except  for  the  6  and  6  incident  angles  that  define  the 
transmitter  location  with  respect  to  the  object. 


Figure  53.  Calculate  Bistatic  RCS  GUI 


For  an  incidence  angle  of  45  degrees  in  the  $  =  0°  plane  and  the  observation 
direction  varied  from  0  to  360  degrees  in  the  <f>  =  0°  plane,  using  a  frequency  of  0.2  GHz 
and  the  number  of  Taylor  Series  terms  set  at  3,  the  bistatic  RCS  of  the  box  model  is 
shown  in  Figure  54.  The  lobe  at  135  degrees  is  the  specular  reflection  from  the  top.  The 

forward  scatter  is  also  evident.  Note  that  the  forward  scatter  is  higher  than  the  highest 
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monostatic  backscatter  lobe.  The  accompanying  ray  diagram  in  Figure  53  illustrates  the 
directions  of  the  expected  large  RCS  lobes. 


Status  GUIs  similar  to  that  of  the  Calculate  Monostatic  RCS  module  were 
designed  to  enable  the  user  to  easily  correct  specific  parameters.  The  Calculate  Bistatic 
RCS,  however  has  an  additional  GUI  for  the  incidence  angle  as  shown  in  Figure  56. 
Other  helpful  tips  may  be  found  in  Calculate  Bistatic  RCS  Help  GUI  shown  in  Figure  57. 


#  Incident  Angle  Status 


m 


© 


Enter  an  incident  angle  between  -360  and  +360  degrees. 


OK 


Figure  56.  Incident  Angle  Status  GUI 


BISTATIC  RCS  Help 


Calculate  Bistatic  RCS  computes  the  total  radar  cross  section  of  the  target  P 


using  a  bistatic  radar  (transmitter  and  receiver  at  different  locations). 

Incidence  Angle  -  The  incidence  direction  is  fixed  and  is  defined  by  a  theta 
and  phi  angle  in  degrees. 

Observation  Angles  -  "WHe  the  incidence  angle  is  fixed,  the  observation 
direction  changes.  Theta/Phi  Start  and  Stop  Angles  are  measured  in  degrees 
and  entering  a  phi  start  and  phi  stop  with  equal  values  results  in  a  phi-cut. 
Similarly,  entering  a  theta  start  and  theta  stop  with  equal  values  results  in  a 
theta-cut. 

Increment  -  defines  the  number  of  degrees  (or  fraction  of  a  degree)  increase 
in  step  for  the  thetafchi  start  and  stop  angles  to  be  used  in  calculating  the 
bistatic  radar  cross  section  of  the  model. 

Surface  Roughness  -  A  model  may  be  designed  to  have  a  rough  surface  by 
inputting  values  for  correlation  distance  and  standard  deviatioa  These 
parameters  contribute  a  disuse  scattering  component  of  the  RCS  and  is  adde 


Close 


Figure  57.  Calculate  Bistatic  RCS  Help  GUI 
59 


THIS  PAGE  INTENTIONALLY  LEFT  BLANK 


60 


V.  SUMMARY  AND  RECOMMENDATIONS 


A.  SUMMARY 

The  main  purpose  of  this  work  was  to  develop  a  Graphical  User  Interface  (GUI) 
that  will  provide  a  convenient  way  of  modeling  complex  objects  using  triangular  facets, 
provide  options  to  define  the  surface  roughness  of  the  scattering  object,  and  calculate  the 
monostatic  or  bistatic  radar  cross  section  over  a  range  of  frequencies  or  wavelengths  of 
interest.  Care  was  taken  to  design  the  GUI  to  make  it  appealing  and  easy  to  use  and 
several  validation  techniques  were  incorporated  to  make  it  robust  and  prevent  the  user 
from  encoding  erroneous  entries.  Help  windows  and  messages  guide  the  user  in  making 
efficient  use  of  the  GUI.  The  GUI  was  implemented  in  MATLAB  utilizing  its  scientific 
computing  features  and  several  built-in  plotting  functions  to  visualize  the  model 
geometry  and  results  of  the  RCS  calculation.  Implementing  the  code  as  a  user  interface 
is  the  best  way  to  avoid  errors  that  may  have  significant  consequences  related  to  the  use 
of  this  tool  for  the  desired  purpose. 

The  GUI  can  be  used  by  both  students  and  electromagnetic  professionals  to 
quickly  determine  the  RCS  performance  of  their  model  by  varying  the  computational 
parameters  to  suit  their  design  requirements. 

B.  RECOMMENDATIONS 

The  GUI  developed  was  based  on  the  basic  physical  optics  approximation.  An 

important  area  for  future  work  is  to  be  able  to  upgrade  this  GUI  into  one  that  incorporates 

other  scattering  mechanisms,  primarily  the  effects  of  multiple  reflections.  Multiple 

reflections  complicate  the  PO  approximation  because  the  incident  field  on  a  patch  is  not 

only  determined  by  the  direct  field  from  the  source,  but  also  the  field  scattered  by  other 
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facets  in  the  direction  of  the  triangle  under  consideration.  Numerical  integrations  will  be 
required  that  substantially  increase  the  computation  time. 
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APPENDIX  A.  MATLAB  M-FILES 


The  following  are  the  description  of  the  Matlab  script  files  contained  in  the 
POFACETS  Graphical  User  Interface: 

1.  BPlot.m 

This  file  generates  either  a  theta  of  phi-cut  bistatic  RCS  calculation  plot  of  a 
selected  model. 

2.  BPlotl.m 

This  file  generates  a  contour  plot  of  the  bistatic  RCS  calculation  of  a  selected 

model. 

3.  BMessage.m 

This  file  displays  the  name  of  the  model  selected  in  the  Calculate  Bistatic  RCS 
GUI  and  enables  the  Calculate  RCS  button. 

4.  G.m 

This  file  calculates  the  recursive  Green’s  function  using  the  two  parameters  (n  and 
y)  passed. 

5.  MPIot.m 

This  file  generates  either  a  theta  of  phi-cut  monostatic  RCS  calculation  plot  of  a 
selected  model. 

6.  MPlotl.m 

This  file  generates  a  contour  plot  of  the  monostatic  RCS  calculation  of  a  selected 

model. 
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7. 


MaxPoints.m 


This  file  displays  a  message  that  the  GUI  design  can  accommodate  only  a 
maximum  of  40  vertices  only. 

8.  Repl_File.m 

This  file  displays  a  GUI  prompting  the  user  to  replace  the  data  of  an  existing 

model. 

9.  Save_As.m 

This  file  displays  a  GUI  with  a  list  of  the  existing  models  and  prompts  the  user  to 
enter  a  filename. 

10.  SvAs.m 

This  file  implements  the  functionality  of  the  Save_As  GUI  by  displaying  the 
filename  selected  from  the  list  and/or  checking  if  the  filename  entered  already  exists. 

11.  VertS  tatus.m 

This  file  displays  a  GUI  that  prompts  the  user  to  check  the  number  of  vertices  if 
the  value  entered  is  erroneous. 

12.  bistatic.m 

This  file  displays  the  Calculate  Bistatic  RCS  GUI. 

13.  bistat.m 

This  file  implements  the  functionalities  of  the  Calculate  Bistatic  RCS  GUI  by 
checking  all  the  parameters  entered,  displaying  GUIs  with  appropriate  messages, 
calculating  the  bistatic  RCS  of  the  selected  model,  and  generating  the  plot. 
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14.  design.m 

This  file  displays  the  Design  Model  Coordinates  GUI  for  encoding  the 
coordinates  of  the  model  being  designed  or  edited. 

15.  dsgn.tn 

This  file  implements  the  functionalities  of  the  Design  Model  Coordinates  GUI  by 
checking  the  individual  coordinates  entered  and  the  whole  coordinate  set. 

16.  dsgnFN.m 

This  file  clears  previous  design  GUIs  and  displays  a  new  Design  Model 
Coordinates  GUI. 

17.  dsgnFO.m 

This  file  displays  the  Open  Model  File  GUI  showing  the  list  of  existing  models 
which  the  user  may  select  to  display  the  data  set  and  plot  of  the  geometry. 

18.  facets.m 

This  file  displays  the  Design  Model  Facets  GUI  for  defining  the  configuration  of 
the  component  triangular  facets  of  the  model. 

19.  fact.m 

This  file  calculates  the  factorial  of  a  number  n  given  the  parameter. 

20.  fcets.m 

This  file  implements  the  functionalities  of  the  Design  Model  Facets  GUI  by 
checking  the  number  of  facets,  individual  nodes  and  set  of  nodes,  illumination  and 
resistivity  values. 
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21.  hlpbistattxt.m 

This  file  contains  the  text  that  is  displayed  in  the  Calculate  Bistatic  RCS  Help 

GUI. 

22.  hlpcoordtxt.m 

This  file  contains  the  text  that  is  displayed  in  the  Design  Model  Coordinates  Help 

GUI. 

23.  hlpfacettxt.m 

This  file  contains  the  text  that  is  displayed  in  the  Design  Model  Facets  Help  GUI. 

24.  hlpmaintxt.m 

This  file  contains  the  text  that  is  displayed  in  the  POFACETS  Main  GUI. 

25.  hlpmodeltxt.m 

This  file  contains  the  text  that  is  displayed  in  the  Design  Model  Geometry  Help 

GUI. 

26.  hlpmonotxt.m 

This  file  contains  the  text  that  is  displayed  in  the  Calculate  Monostatic  RCS  Help 

GUI. 

27.  hlpbistabm 

This  file  displays  the  Calculate  Bistatic  RCS  Help  GUI. 

28.  hlpdsgncoord.m 

This  file  displays  the  Design  Model  Coordinates  Help  GUI. 

29.  hlpdsgnfacet.m 

This  file  displays  the  Design  Model  Facets  Help  GUI. 


66 


30.  hlpmain.m 

This  file  displays  the  POFACETS  Main  Help  GUI. 

31.  hlpmodel.m 

This  file  contains  the  text  that  is  displayed  in  the  Design  Model  Facets  Help  GUI. 

32.  hlpmono.m 

This  file  displays  the  Calculate  Monostatic  RCS  Help  GUI. 

33.  message.m 

This  file  displays  the  name  of  the  model  selected  in  the  Calculate  Monostatic  RCS 
GUI  and  enables  the  Calculate  RCS  button. 

34.  mod.m 

This  file  implements  the  functionalities  of  the  Model  Geometry  GUI  by  plotting 
the  model,  labeling  its  vertices  and  facets,  displaying  the  grid,  and  printing  the  data  set 
and  model  plot. 

35.  model.m 

This  file  displays  the  Model  Geometry  GUI. 

36.  mono.m 

This  file  implements  the  functionalities  of  the  Calculate  Monostatic  RCS  GUI  by 
checking  all  the  parameters  entered,  displaying  GUIs  with  appropriate  messages, 
calculating  the  monostatic  RCS  of  the  selected  model,  and  generating  the  plot. 

37.  monostatic.m 

This  file  displays  the  Calculate  Monostatic  RCS  GUI. 


67 


r 


38.  mplt.m 

This  file  prints  the  Monostatic  or  Bistatic  RCS  plot. 

39.  msgComputing.m 

This  file  displays  a  message  that  the  GUI  is  calculating  either  the  monostatic  or 
bistatic  RCS  of  the  model  selected. 

40.  pofacets.m 

This  file  displays  the  POFACETS  Main  GUI  where  the  user  may  choose  to  run 
any  of  the  three  modules  -  Design  Model,  Calculate  Monostatic  RCS  and  Calculate 
Bistatic  RCS. 
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